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Abstract 



We calculate the pressure for pure-glue QCD at high temperature to two-loop order 
using hard-thermal-loop (HTL) perturbation theory. At this order, all the ultraviolet 
divergences can be absorbed into renormalizations of the vacuum energy density and 
the HTL mass parameter. We determine the HTL mass parameter by a variational 
prescription. The resulting predictions for the pressure fail to agree with results from 
lattice gauge theory at temperatures for which they are available. 
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I. INTRODUCTION 



Relativistic heavy-ion collisions allow the experimental study of hadronic matter at 
energy densities exceeding that required to create a quark-gluon plasma. A quantitative 
understanding of the properties of a quark-gluon plasma is essential in order to determine 
whether it has been created. Because QCD is asymptotically free, its running coupling 
constant a s becomes weaker as the temperature increases. One might therefore expect 
the behavior of hadronic matter at sufficiently high temperature to be calculable using 
perturbative methods. Unfortunately, a straightforward perturbative expansion in powers 
of a s does not seem to be of any quantitative use even at temperatures orders of magnitude 
higher than those achievable in heavy-ion collisions. 

The problem is evident in the free energy T of the quark-gluon plasma, whose weak- 
coupling expansion has been calculated through order 0PJ 2 fT)-^. For a pure-glue plasma, 
the first few terms in the weak-coupling expansion are 



J~ QCD — ^ideal 



15 a s / a x \ 3 / 2 
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(i) 



where J-ideai = — (87r 2 /45)T 4 is the free energy of an ideal gas of massless gluons and 
a s = a s (fi) is the running coupling constant in the MS scheme. In Fig. [|, the free en- 
ergy is shown as a function of the temperature T/T c , where T c is the critical temperature 
for the deconfinement transition. The weak-coupling expansions through orders a s , a^ 2 , 
a 2 , and oPj 2 are shown as bands that correspond to varying the renormalization scale \i by 
a factor of two from the central value /i = 2ttT. As successive terms in the weak-coupling 
expansion are added, the predictions change wildly and the sensitivity to the renormaliza- 
tion scale grows. It is clear that a reorganization of the perturbation series is essential if 
perturbative calculations are to be of any quantitative use at temperatures accessible in 
heavy-ion collisions. 

The free energy can also be calculated nonperturbatively using lattice gauge theory [|J . 
The thermodynamic functions for pure-glue QCD have been calculated with high precision 
by Boyd et al. ||. There have also been calculations with Nf = 2 and 4 flavors of dynamical 
quarks ||. In Fig. [I], the lattice results for the free energy of pure-glue QCD from Boyd et 
al. |5j are shown as diamonds. The free energy is very close to zero near T c . As the tem- 
perature increases, the free energy increases and approaches that of an ideal gas of massless 
gluons. We will regard the lattice results as the correct results for the thermodynamic func- 
tions. One goal of any reorganization of perturbation theory is to obtain a free energy that 
agrees within its domain of validity with the lattice results. 

There is of course little to be gained by just reproducing the results of lattice gauge 
theory. A method for reorganizing perturbation theory is of practical use only if it allows the 
calculation of quantities that are not so easily calculated using lattice gauge theory. There 
are many observables that are difficult or even impossible to calculate using lattice gauge 
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T/T c 

FIG. 1. The free energy for pure-glue QCD as a function of T/T c . The weak-coupling expan- 
sions through orders a s , c?J 2 , a 2 s , and a J 2 are shown as bands that correspond to varying the 
renormalization scale /i by a factor of two. The diamonds are the lattice result from Boyd et al. || . 
The size of the diamonds indicate the typical error bar. 



theory. First, lattice gauge theory becomes increasingly inefficient at higher temperatures, 
so some other method is required to extrapolate to high T. Second, calculations with light 
dynamical quarks require orders of magnitude more computer power than pure-glue QCD. 
Third, the Monte Carlo approach used in lattice gauge theory fails completely at nonzero 
baryon number density. Finally, lattice gauge theory is only effective for calculating static 
quantities, but many of the more promising signatures for a quark-gluon plasma involve 
dynamical quantities. 

The only rigorous method available for reorganizing perturbation theory in thermal QCD 
is dimensional reduction to an effective 3-dimensional field theory [/],§]. The coefficients of 
the terms in the effective lagrangian are calculated using perturbation theory, but calcula- 
tions within the effective field theory are carried out nonperturbatively using lattice gauge 
theory. Dimensional reduction has the same limitations as ordinary lattice gauge theory: 
it can be applied only to static quantities and only at zero baryon number density. Un- 
like in ordinary lattice gauge theory, light dynamical quarks do not require any additional 
computer power, because they only enter through the perturbatively calculated coefficients 
in the effective lagrangian. This method has been applied to the Debye screening mass for 
QCD as well as the pressure 0. 

There are some proposals for reorganizing perturbation theory in QCD that are es- 
sentially just mathematical manipulations of the weak coupling expansion. The methods 
include Pade approximates ||, Borel resummation [ID||, and self-similar approximates [ |TT 



These methods are used to construct more stable sequences of successive approximations 
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that agree with the weak-coupling expansion when expanded in powers of a s . These methods 
can only be applied to quantities for which several orders in the weak-coupling expansion 
are known, so they are limited in practice to the thermodynamic functions. 

One promising approach for reorganizing perturbation theory in thermal QCD is to use 
a variational framework. The free energy T is expressed as the variational minimum of a 



thermodynamic potential fi(T, a s \ 
that we denote collectively by m 2 : 



m 



that depends on one or more variational parameters 



T{T, a s ) = Q(T,a s ;m 2 



dn/dm 2 =0 



(2) 



A particularly compelling variational formulation is the $- derivable approximation, in which 
the complete propagator is used as an infinite set of variational parameters The $- 

derivable thermodynamic potential Q is the 2PI effective action, the sum of all diagrams 
that are 2-particle-irreducible with respect to the complete propagator [DJ. The n-loop 
^-derivable approximations, in which Q is the the sum of 2PI diagrams with up to n loops, 
form a systematically improvable sequence of variational approximations. Until recently, $- 
derivable approximations have proved to be intractable for relativistic field theories except 
for simple cases in which the self-energy is momentum- independent. However there has 
been some recent progress in solving the 3-loop ^-derivable approximation for scalar field 
theories. Braaten and Petitgirard have developed an analytic method for solving the 3- 

6 4 field theory @ 



loop ^-derivable approximation for the massless 



Van Hees and Knoll 



have developed numerical methods for solving the 3-loop ^-derivable approximation for the 
b A field theory [ITH 



massive 



They have also investigated renormalization issues associated 
with the ^-derivable approximation. 

The application of the ^-derivable approximation to QCD was first discussed by McLer- 
ran and Freedman |16| . One problem with this approach is that the thermodynamic potential 
Q is gauge dependent, and so are the resulting thermodynamic functions. The gauge depen- 
dence is the same order in a s as the truncation error. However the most serious problem is 
that even the 2-loop ^-derivable approximation has proved to be intractable. 

The 2-loop ^-derivable approximation for QCD has been used as the starting point for 



HTL resummations of the entropy by Blaizot, Iancu and Rebhan |17| and of the pressure 
by Peshier []18| . The thermodynamic potential ^2-loop is a functional of the complete gluon 
propagator D^P). The HTL resummations of Refs. [17| and [IS] can be derived in 2 
steps. The first step is to replace the 2-loop functional at its variational point by a 1-loop 
functional evaluated at the 2-loop variational point. In the resummation of the pressure 



of Ref . 



the 2-loop functional is the thermodynamic potential and this step is a weak- 



coupling approximation: 



n 



2— loop 



^l-loop [D i 



2 — loop ~ 



(3) 



2 — loop ~ 



In the resummation of the entropy of Ref. ||17|| , the 2-loop functional is the derivative of 
^2-ioop with respect to T and this step is an exact equality. The second step exploits the 
fact that the HTL gluon propagator D^J L (P) is an approximate solution to the variational 



equation 5fi 2 - 



loop 



0. The HTL gluon propagator depends on one parameter m D , which 
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can be interpreted as the Debye screening mass for the gluon. The HTL gluon propagator 
satisfies the variational equation to leading order in a s provided that m 2 D reduces in the 
weak-coupling limit to 

2 ^N c 2 

m D = — — a.(ji)T , (4) 

with some appropriate choice for the scale /i such as fi = 2ttT . Thus we can approximation 
the solution to the variational equation in (j3|) by D^J h (P): 



HTL 1 



^l-loop[-D/ii/] ~ ^l-loop[-D M ^ 

5n2_ioop=o 



(5) 

=4ttq s T 2 



This approximate solution holds when m 2 D is given by (Q), however, there is some freedom 
in the choice of the parameter m 2 D , as long as it reduces to (f|) in the weak-coupling limit. 
It can not be determined variationally, because the variational character of the thermody- 
namic potential was lost in the first step ([3]). With the prescription (f|), the errors in the 
thermodynamic functions are of order o?J 2 . The errors can be reduced to order a 2 by adding 
an o?J 2 term to the right side of (f|). 

The intractability of ^-derivable approximations motivates the use of simpler variational 
approximations. One such strategy that involves a single variational parameter m has been 
called optimized perturbation theory , variational perturbation theory ]20| , or the linear 



5 expansion [^T|. This strategy was applied to the thermodynamics of the massless 4 
field theory by Karsch, Patkos and Petreczky under the name screened perturbation theory 
|22fl . The method has also been applied to spontaneously broken field theories at finite 



temperature |[23|| . The calculations of the thermodynamics of the massless 4 field theory 



using screened perturbation theory have been extended to 3 loops |24| . The calculations can 



be greatly simplified by using a double expansion in powers of the coupling constant and 
m/T 

HTL perturbation theory (HTLpt) is an adaptation of this strategy to thermal QCD [EB 



The exactly solvable theory used as the starting point is one whose propagators are the HTL 
gluon propagators. The variational mass parameter can be identified with the Debye 
screening mass. The one-loop free energy in HTLpt was calculated for pure-glue QCD in 
Ref. [26f and for QCD with light quarks in Ref. ||27|| . At this order, the parameter mjj cannot 



be determined variationally, so the prescription (H) was used. The resulting thermodynamic 
functions have errors of order a s , but the terms of order cfij 2 associated with Debye screening 
are correct. A two- loop calculation is required to reduce the errors to order a 2 . At two-loop 
order, it is also possible to determine using a variational prescription. 



One difference between HTLpt and the HTL resummation methods of Refs. [|T7| and [fTS] 
is in how they deal with gauge invariance. HTLpt is constructed in such a way that physical 
observables are gauge invariant order-by-order in perturbation theory. Gauge invariance 
arises in the same way as in ordinary perturbation theory by cancellations between diagrams. 
In the HTL resummation methods of Refs. |L7| and [|l8j, the 2- loop thermodynamic potential 



^2-ioop that is used as the starting point is gauge dependent. In the first step ([3]) of the 
derivation, f2 2 _i oop is replaced by a 1-loop functional f2i_i oop that is gauge invariant, but 
the variational equation £fi2-ioop = is still gauge dependent. In the second step the 
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solution to that variational equation is approximated by D^J L , and it is only at this point 



that the gauge dependence disappears. 

Another difference between HTLpt and the HTL resummation methods of Refs. |I7J and 
18fl is in the ranges of observables to which they can be applied. The HTL resummation 



methods were specifically formulated as approximations to the thermodynamic functions, so 
they cannot be easily applied to other observables. However, they can be used to calculate 
the thermodynamic functions in cases where calculations using conventional lattice gauge 
theory are difficult or impossible: the high-temperature limit of pure-glue QCD, QCD with 
light quarks, and QCD with nonzero baryon number density. In contrast to these methods, 
HTLpt has the same wide range of applicability as ordinary perturbation theory. It can be 
used to calculate the thermodynamic functions, but it can also be applied to all the standard 
signatures of a quark-gluon plasma. It has some of the limitations of ordinary perturbation 
theory. Calculations can be carried out only up to the order at which the magnetic screening 
problem causes diagrammatic methods to break down. 

In this paper, we calculate the thermodynamic functions of QCD to 2-loop order in 
HTLpt. We begin with a brief summary of HTLpt in Section [D. in Section [TTT| , we give 
the expressions for the one-loop and two-loop diagrams for the thermodynamic potential. 
In Section [TV], we reduce those diagrams to scalar sum-integrals. We are unable to compute 
those sum-integrals, so in Section |V|, we evaluate them approximately by expanding them 
in powers of rno/T. The diagrams are combined in Section |Vl]to obtain the final results for 
the two- loop thermodynamic potential up to 5 th order in g and m D /T. In Section |V1J| , we 
present our numerical results for the thermodynamic functions of QCD. There are several 
appendices that contain technical details of the calculations. In Appendix we give the 
Feynman rules for HTLpt in Minkowski space to facilitate the application of this formalism 
to signatures of the quark-gluon plasma. The most difficult aspect of these calculations 
was the evaluation of the sum- integrals obtained from the expansion in mn/T. We give 
the results for these sum-integrals in Appendix |B|. The evaluation of some difficult thermal 
integrals that were required to obtain the sum-integrals is described in Appendix |C[ 



II. HTL PERTURBATION THEORY 

The lagrangian density that generates the perturbative expansion for pure-glue QCD 
can be expressed in the form 

£qcd = ~ -Tr (G^ U G^ U ) + £ g f + £ g host + A£q CD , (6) 

where G^ v = d fM A u — d v A^ — ig[A^ A v ] is the gluon field strength and A^ is the gluon 
field expressed as a matrix in the SU(N C ) algebra. The ghost term £ g host depends on the 
choice of the gauge-fixing term £ g f. Two choices for the gauge- fixing term that depend on 
an arbitrary gauge parameter £ are the general covariant gauge and the general Coulomb 
gauge: 

£ gf = -^Tr ((<9^) 2 ) covariant , (7) 

= --Tr ((V • A) 2 ) Coulomb . (8) 
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The perturbative expansion in powers of g generates ultraviolet divergences. The renormal- 
izability of perturbative QCD guarantees that all divergences in physical quantities can be 
removed by renormalization of the coupling constant a s = g 2 /47r. If we use dimensional reg- 
ularization with minimal subtraction as a renormalization prescription, the renormalization 
can be accomplished by substituting a s — > a s + Aa s , where the counterterm Aa s is a power 
series in a s whose coefficients have only poles in e: 

11JV C o / 121iV c 2 17iV 2 \ 3 4N 
127re \1447r z e z 48n z e ) 

Renormalized perturbation theory can be implemented by including among the interactions 
terms a counterterm lagrangian A£qcd that is given by the change in the first 3 terms on 
the right side of (H) upon substituting g — > g(l + Ao^) 1 / 2 . 

Hard-thermal-loop perturbation theory (HTLpt) is a reorganization of the perturbation 
series for thermal QCD. The lagrangian density is written as 

C = (£qcd + £ HT l) \ g ^ g + A£ H tl. (10) 
The HTL improvement term is 



Am. = ~ S)mlTt (g/mx [jf^y ) G ^ I » ( n ) 



where is the covariant derivative in the adjoint representation, = (1, y) is a light-like 
four- vector, and (. . .) y represents the average over the directions of y. The term flTT| ) has the 
form of the effective lagrangian that would be induced by a rotationally invariant ensemble 
of colored sources with infinitely high momentum. The parameter can be identified with 
the Debye screening mass. HTLpt is defined by treating 5 as a formal expansion parameter. 
The free lagrangian in general covariant gauge is obtained by setting 5 = in (|10D : 



£ free = -Tr (d^A^A" - dpAvVA") - ^Tr ((^) 



-im^Tr Ud,A a - d a A,) ) {d»A p - d p A») } . (12) 



' v 



The resulting propagator is the HTL gluon propagator. The remaining terms in (0) are 
treated as perturbations. The Feynman rules for gluon and ghost propagators and the 
3-gluon, ghost-gluon, and 4-gluon vertices are given in Appendix 0. 

The HTL perturbation expansion generates ultraviolet divergences. In QCD perturba- 
tion theory, renormalizability constrains the ultraviolet divergences to have a form that can 
be cancelled by the counterterm lagrangian A£qcd- There is no proof that the HTL per- 
turbation expansion is renormalizable, so the general structure of the ultraviolet divergences 
is not known. The most optimistic possibility is that HTLpt is renormalizable, so that the 
ultraviolet divergences in physical quantities can all be cancelled by renormalization of the 
coupling constant a s , the mass parameter m 2 D , and the vacuum energy density Sq. If this 
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is the case, the renormalization of a physical quantity can be accomplished by substituting 
a s — > a s + Aa s and m 2 D — > m 2 D + Am 2 D , where Aa s and Am 2 D are counterterms. In the case 
of the free energy, it is also necessary to add a vacuum energy counterterm A£ - If we use 
dimensional regularization with minimal subtraction as a renormalization prescription, the 
form of the counterterms for 5a s , (1 — 5)m D , and Sq should be the power of (1 — S)m 2 D re- 
quired by dimensional analysis multiplied by a power series in 5a s with coefficients that have 
only poles in e. The counterterm for Sa s should be identical to that in ordinary perturbative 
QCD given in @ with: 

11N c, 2 2 f 121N 2 17N 2 \ ^ 3 _ 4 . 
SAa ° = -12^ a ' + 145* - 48^ SV - + °^ ' (13) 



The leading term in the delta expansion of the £q counterterm A£ was deduced in Ref. [26 
by calculating the free energy to leading order in 5. The So counterterm A£ must therefore 
have the form 

/ 7V 2 - 1 \ 
A£o = + 0(5a s )j (1 - 5) 2 m D . (14) 

To calculate the free energy to next-to-leading order in 5, we need the counterterm ASq 
to order 5 and the counterterm Am 2 D to order 8. We will show that there is a nontrivial 
cancellation of the ultraviolet divergences if the mass counterterm has the form 

Ami = Sa s + 0{5 2 a*)) (1 - 5)m 2 D . (15) 

Renormalized perturbation theory can be implemented by including a counterterm la- 
grangian A£htl among the interaction terms in (|10|). 

Physical observables are calculated in HTLpt by expanding them in powers of 5, trun- 
cating at some specified order, and then setting 5=1. This defines a reorganization of the 
perturbation series in which the effects of the m 2 D term in (|TJ) are included to all orders 
but then systematically subtracted out at higher orders in perturbation theory by the 8m 2 D 
term in (0). If we set 5 = 1, the lagrangian ( |T0l) reduces to the QCD lagrangian (|). If 
the expansion in 5 could be calculated to all orders, all dependence oiiuid should disappear 
when we set 5 = 1. However, any truncation of the expansion in 5 produces results that 
depend on mo- Some prescription is required to determine mo as a function of T and a s . 
We choose to treat as a variational parameter that should be determined by minimizing 
the free energy If we denote the free energy truncated at some order in 5 by Q(T, a s , mo, 8), 
our prescription is 

° n(T,a s ,m D ,6= 1) = 0. (16) 



dm o 

Since fi(T, a s , mo, 5 = 1) is a function of a variational parameter mo, we will refer to it 
as the thermodynamic potential. We will refer to the variational equation (|T6|) as the gap 
equation. The free energy T is obtained by evaluating the thermodynamic potential at the 
solution to the gap equation. Other thermodynamic functions can then be obtained by 
taking appropriate derivatives of JF with respect to T . 
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III. DIAGRAMS FOR THE THERMODYNAMIC POTENTIAL 



The thermodynamic potential at leading order in HTL perturbation theory (HTLpt) is 



^lo = (N 2 C - \)T 9 + A £ t 



o • 



(17) 



where T g is the contribution to the free energy from each of the color states of the gluon: 



.F s = -^{(d-l)log[-A T (P)] + logA L (P)} . 



(18) 



The transverse and longitudinal HTL propagators A^(P) and A^(P) are given in (|A.49|) 
and (|A.50|) . We use dimensional regularization with d = 3 — 2e spatial dimensions to 
regularize ultraviolet divergences. The term of order 5° in the vacuum energy counterterm 



was determined in Ref. 26 



An£ ( 



o<-o 



1287r 2 e 



m 



D 



(19) 





J 7 

<J A, 



4g 



3g 



•J rt 



J gh 



FIG. 2. Diagrams contributing through NLO in HTLpt. Shaded circles indicate dressed HTL 
propagators and vertices. 



The thermodynamic potential at next-to-leading order in HTLpt can be written 



NLO 



^lo + (iV c 2 - 1) [Fz g + F ig + T gh + J^htl] + A^o + A x m 2 D 



2 9 n 



(20) 



where A\Sq and Aim 2 D are the terms of order 5 in the vacuum energy density and mass 
counterterms. The contributions from the two-loop diagrams with the three-gluon and four- 
gluon vertices are 



N r 
12 Tpq 



N r 

F± g = -^g 2 ± r^< ACT (P, -P, Q, -Q)A^(P)A X °(Q) . 
8 Tpq 



(21) 
(22) 



Expressions for the gluon propagator tensor A Ml/ , the 3-gluon vertex tensor T AtAp , and the 
4-gluon vertex tensor Y tlu,Xcr in Minkowski space are given in ( |A.25[ ) or ( |A.26| ), ( |A.32j ), and 
( |A.41| ). Prescriptions for translating them into the Euclidean tensors appropriate for the 



imaginary time formalism are given in Appendix |A 7|. The contribution from the ghost 
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diagram depends on the choice of gauge. The expressions in the covariant and Coulomb 
gauges are 

T 9h = Y 92 f PQ -^2^ RV ^^ P ) covariant , (23) 

= —9 2 i 44 (Q" -Q-nn»)(R u - R-nn v )W(P) Coulomb. (24) 
2 Tpq q 2 r 2 

The contribution from the HTL counterterm diagram is 

^HTL = ^ p n^(P)A^(P). (25) 

It can also be obtained by substituting m 2 D — > (1 — 5)m 2 D in the one-loop expression T g in 
(|18|) and expanding to first order in 5: 

^htl = \f p [{d - 1)II T (P)A T (P) - n L (P) A L (P)] . (26) 

Provided that HTLpt is renormalizable, the ultraviolet divergences at any order in 6 can 
be cancelled by renormalizations of the vacuum energy density Sq, the HTL mass parameter 
mjj, and the coupling constant a s . Renormalization of the coupling constant does not enter 
until order 5 2 . We will calculate the thermodynamic potential as a double expansion in 
powers of g and ttid/T, including all terms through 5 th order. The 5a s term in ASo does not 
contribute until 6 th order in this expansion, so the term of order 5 in ASq can be obtained 
simply by expanding (TTO) to first order in 5: 



N 2 - 1 „ 

Axfc = --t^n • (27) 

The remaining ultraviolet divergences must be removed by renormalization of the mass 
parameter mo- We will find that there are ultraviolet divergences in the a s m 2 D T 2 and 
a s m 3 D T 3 terms, and both are removed by the same counterterm Aimfj. This provides 
nontrivial evidence for the renormalizability of HTLpt at this order in 5. 



The sum of the 3-gluon, 4-gluon, and ghost contributions in (|2l|) , (|22D , and (|23|) or ( |24|) 



is gauge- invariant. By inserting the expression ( |A.25| ) or (|A.26| ) for the gluon propagator 



tensor and using the Ward identities (|A.35|) and ( |A.42|) , one can easily verify that the sum 



of these three diagrams is independent of the gauge parameter £ in both covariant gauge 
and Coulomb gauge. With more effort, we can verify the equivalence of the covariant gauge 
expression with £ = (Landau gauge) and the Coulomb gauge expression with £ = 0. This 
involves expanding the tensor npHp in the covariant gauge propagator into the sum of terms 
proportional to n fi n u , P fJ, n v n ll P u , and P^P U , and then applying the Ward identities to the 
terms involving P M or P v . 

IV. REDUCTION TO SCALAR SUM-INTEGRALS 

The first step in calculating the thermodynamic potential is to reduce the sum of the 
diagrams to scalar sum-integrals. The one-loop diagram in (111) and the HTL counterterm 
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diagram fl25|) are already expressed in terms of scalar integrals. We proceed to consider the 
3-gluon diagram in (pll), the 4-gluon diagram in (p2l), and the ghost diagram in Landau 
gauge which is given in (^3|). The expression for the sum of these three diagrams is simpler 
than that of the 3-gluon diagram alone. We insert the gluon propagator in the form (|A.29|) 
with £ = 0. It has terms proportional to At and Ax, where Ax is the combination of 
transverse and longitudinal propagators defined in (|A.27 ). When a momentum P M from the 



gluon propagator tensor is contracted with a 3-gluon or 4-gluon vertex, the Ward identities 
can be used to reduce it ultimately to expressions involving the inverse propagator ( |A.20| ). 



The term Ay/A^ can be eliminated in favor of A x /Al using the definition ( |A.27| ). This 
reduces the sum of the 3-gluon, 4-gluon, and ghost diagrams to the following form: 

J~3g + J~&g + J-gh 



—g 2 t jr^ A r^ A A T (P)A T (Q)A T (R) - 3r^°r^° a t (p)a t (Q)a x (r) 

12 Tpq [ 

+3r M>o r M>o a t (P)A x (Q)A x (P) - (r 000 ) 2 A X (P)A X (Q)A X (R) 
+3d(d + 1) A T (P) A T (Q) - 6d A T (P)A X (Q) + h 00 ' 00 A X (P)A X (Q) 

+6 (9± A T (P)- n - Qn - R A x (P)) 
\Q 2 R 2 n ; Q 2 R 2 xy 7 

10 ( n-Q riQ-R n-Q n-R \ A X (Q) 

( n-Qn-Rn Q -n R n-Qn-R \ A X (Q) A X (R) \ 

+6 ( fr* At(P) "^^ Ax(P) J A L (Q) A L (R)j- ( 8) 

In the 3-gluon and 4-gluon vertex tensors, we have suppressed the momentum arguments: 
r ^A = r ^A( Pj g ; P ) and r oo,oo = r oo,oo( P; _ P) q ; _Qy 

The next step is to insert the Euclidean analogs of the expressions (|A.32|) and ( A. 41 ) 



for the 3-gluon and 4-gluon vertex tensors. The combinations of terms that appear in ( |28|) 
can be simplified using the 'Ward identities" (|7Q| ), ( [Oil) , and flA~3^ ) satisfied by the HTL 
correction tensors: 

p^Ap^A = 3(J (p2 + Q2 + R 2j + ^ T ^X T „uX ? ^ 

P^op^o =p 2 + q 2 + 4r 2 + 2d(n-P) 2 + 2d(n-Q) 2 - d(n-R) 2 

+2m 2 D (2T R -T P - T Q ) + m% T^°T^° , (30) 
P^oopMOO = 2q 2 + 2r 2 _ p 2 

-2m 2 D [2T P — Tq — T R + n - (Q — R) T 000 ] + m% T^ 00 T^ 00 , (31) 
(r 000 ) 2 = m 4 D (T 000 ) 2 , (32) 
r 00 ' 00 = -m 2 D T 0000 . (33) 

In the 3-gluon and 4-gluon HTL correction tensors, we have suppressed the momentum 
arguments: T^ A = T^ A (P,Q,P) and T 0000 = T 0000 (P,Q, -P, -Q). We have also used 
the short-hand Tp = — T 00 (P, — P) for the 2-gluon HTL correction tensor. Inserting the 



expressions (|29|)-(|3~3D into (|28|) and eliminating 1/A^(P) in favor of Tp, the reduction to 



scalar integrals is 
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S~3 g + T± a + T t 



~gh 



N, 



PQ 



1 



3dR 2 + -m A D T^T^ 
3 



+ 



+ 



A T (P)A T (Q)A T (R) 

-2q 2 - 4r 2 - 4d(n-Q) 2 + d(n-R) 2 

-Am 2 D (?r - Tq) - mjj T^°T^°] A T (P)A T (Q)A X (R) 
-p 2 + Ar 2 - 2m 2 D n-(Q - R) T 000 

-Am 2 D {T P - T R ) + m A D T^T" 00 ! A T (P)A X (Q)A X (R) 



- l -ml (T 000 ) 2 A X (P)A X (Q)A X (R) 

+d(d + 1) A T (P)A T (Q) - 2d A T (P)A X (Q) 

-\m 2 D T 0000 A X (P)A X (Q) 

+2^A T (P) (l - [q 2 + m 2 D (l - T Q )]A X (Q)) 



2^§^A X (P) (l - [g 2 + m 2 D (l - T Q )]A X (Q)) 
x (l-[ r 2 + m 2 D (l-T R )]A x (R) 
+4^ [g 2 + m 2 D (l - T Q )} A t (P)A x (Q) 
(2^-l)q-r 

- [g + "w- T Q)][ r 



x(l 



r R )]A x {R) 



Q 2 R 2 



4(1 



xA T (P)A x (Q)A x (i?) 



(34) 



V. EXPANSION IN THE MASS PARAMETER 

The thermodynamic potential has been reduced to scalar sum-integrals. In Ref. |2T) 



the sum-integrals for the one-loop free energy were evaluated exactly by replacing the sums 
by contour integrals, extracting the poles in e, and then reducing the momentum integrals 
to integrals that were at most two-dimensional and could therefore be easily evaluated 
numerically It was also shown that the sum-integrals could be expanded in powers oimp/T, 
and that the first few terms in the expansion gave a surprisingly accurate approximation to 
the exact result. 

If we tried to evaluate the 2-loop HTL free energy exactly, there are terms such as those 
involving , ]~ fJ-uX T fJ-uX that could at best be reduced to 5-dimensional integrals that would have 
to be evaluated numerically. We will therefore evaluate the sum-integrals approximately by 
expanding them in powers of mo/T. We will carry out the mp/T expansion to high enough 
order to include all terms through order g 5 if m^/T is taken to be of order g. 
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A. 1-loop sum-integrals 

The one-loop sum-integrals include the leading order free energy given by the sum- 
integrals flTgp and the HTL counterterm given by fl2"S|). The leading order free energy must 
be expanded to order (m^/T) 5 in order to include all terms through order g 5 . The HTL 
counterterm has an explicit factor of m 2 Dl so the sum-integral for the HTL counterterm 
diagram need only to be expanded to order (mpt/T) 3 to include all terms through order g 5 . 

The sum- integrals over P involve two momentum scales: rrirj and T. In order to expand 
them in powers of rrio/T, we separate them into contributions from hard loop momentum, 
for which some of the components of P are of order T, and soft loop momenta, for which all 
the components of P are of order mp. We will denote these regions by (h) and (s). Since 
the Euclidean energy P is an integer multiple of 27rT, the soft region requires Pq = 0. 



1. Hard contributions 



If P is hard, the denominators P 2 + and p 2 + in the propagators are of order 
T, but the self-energy functions and are of order m? D . The m^/T expansion can 
therefore be obtained by expanding in powers of and n^. 

For one- loop free energy, we need to expand to second order in m 2 D : 



d 



4(d- 1 



■m 



D 



2\2 



P 2 P 2 



2d-T P 

^4 



p2 p2 



Tp + d—(T P 

pi 



(35) 



Note that the function Tp cancels from the m 2 D term because of the identity ( |A.12| ). The 
values of the sum-integrals are given in Appendix 0. Inserting those expressions, the hard 
contributions to the leading-order free energy reduce to 



CO 



— T + — 

45 24 



1+ 2 + 2 



c(-i), 



1287T 2 



1 n 2tt 2 " 
7 -7 + 2 7 + - 



\4rrTJ 



(JL) 

2e 



2c 



m 2 D T 2 



711 



D 



(36) 



Note that the pole in the m A D term is cancelled by the counterterm (|i~9|). 

The HTL counterterm diagram has an explicit factor of m 2 Dl so we need only to expand 
the sum-integral to first order in m 2 D . Eliminating Ut(P) and Ul(P) in favor of the function 
T P , the result is 



T {h) 



ct 



— m 



D 



pP 2 



2(d 



-m 



D 



1 



1 



(P 



2\2 



p 2 P 2 



P 



p 2 P 2 



P 



(37) 



The values of the sum-integrals are given in Appendix 0. Inserting those expressions, the 
hard contributions to the HTL counterterm in the free energy reduce to 
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Note that the first term in ( j38l) cancels the order-e° term in the coefficient of m 2 D T 2 in (j36[) . 
We have kept the order-e term in the coefficient of m 2 D T 2 in (pop, because it will contribute 
to the final result through the mass counterterm. 



2. Soft contributions 



The soft contribution comes from the Pq = term in the sum-integral. At soft momen- 
tum P = (0, p), the HTL self-energy functions reduce to Ut(P) = and ITl(P) = m 2 D . The 
transverse term vanishes in dimensional regularization because there is no momentum scale 
in the integral over p. Thus the soft contribution comes from the longitudinal term only. 

The soft contribution to the leading order free energy is 



= \T J log(p 2 + ml) . 



Using the expression for the integral in Appendix [Cj we obtain 



127T 



\2m 



D 



2c 



m D T. 



(39) 



(40) 



The soft contribution to the HTL counterterm is 



J ct 



-m D T 



2 " JppH m z D 
Using the expression for the integral in Appendix O, we obtain 

- J_rn 3 T 



(41) 



(42) 



B. 2-loop sum-integrals 

The sum of the two-loop sum-integrals is given in (|3~3|). Since these integrals have an 
explicit factor of g 2 , we need only expand the sum-integrals to order (m^/T) 3 to include all 
terms through order g 5 . 

The sum-integrals involve two momentum scales: mo and T. In order to expand them 
in powers of mrj/T, we separate them into contributions from hard loop momenta and soft 
loop momenta. This gives 3 separate regions which we will denote (hh), (hs), and (ss). In 
the (hh) region, all 3 momenta P, Q, R are hard. In the (hs) region, two of the 3 momenta 
are hard and the other is soft. In the (ss) region, all 3 momenta are soft. 
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1. Contributions from (hh) region 



If P, Q, R are all hard, we can obtain the m^/T expansion simply by expanding in 
powers of m 2 D . To obtain the expansion through order m|,/T 3 , we need only expand to first 
order in m 2 D , with Ax and Hr taken to be of order m 2 D : 



9 2 f pQ UdR 2 A T (P)A T (Q)A T (R) 



-2q 2 - 4r 2 - Ad(n-Q) 2 + d(n-R) 2 ] A T (P)A T (Q)A X (R) 
+d(d + 1) A T (P)A T (Q) - 2d A T (P)A X (Q) 

+ 2 ^A T (P) (l-2q 2 A x (Q) 



-2 



Q 2 R 2 



(43) 



For hard momenta, the self-energies are suppressed by m 2 D /T 2 relative to the propagators, 
so they can be expanded in powers of Hr and IIl. Expanding all terms to first order in m 2 D , 
and using ( |A.6| ) and ( |A.7| ) to eliminate Ut{P) and ITl(P) in favor of Tp, we obtain 



T {hh) _ N c 2 r 
S 3g+ ,g +gh - A 9^ pQ 



1 1 



Nc 2 2 

+ -7-9 ™ D , 
A JPQ 



si 1 , , si 1 



+2 



P 2 (Q 2 ) 2 
' + (d + 2) 



P 2 q 2 Q 2 



P 2 Q 2 R 2 



-2d P ' Q „, - Ad 



p2Q2 r 2 
2 



Q 



P 2 Q 2 (r 



2\2 



A 



P 2 Q 2 (r 2 ) 2 P 2 Q 2 r 2 R 2 



+Ad- 



P 2 Q 2 (r 2 ) 2 n P 2 Q 2 (r 2 ) 2 
Inserting the sum-integrals from Appendix EL this reduces to 



(44) 



{hh) 



Zg+Ag+gh u 



96 



1 



4.621 



N c a s ( \i 



3tt V4vrT 



Ac 



2 ml 

m D T 



(45) 



2. The (hs) contributions 

In the (hs) region, the soft momentum can be any one of the three momenta P, Q, or 
R. However we can always permute the momenta so that the soft momentum is P = (0, p). 
The function that multiplies the soft propagator A T (0, p) or A x (0, p) can be expanded in 
powers of the soft momentum p. In the case of Ay(0, p), the resulting integrals over p 
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have no scale and therefore vanish in dimensional regularization. The integration measure 
/ scales like m 3 D , the soft propagator Ax(0, p) scales like l/m 2 D , and every power of p in 
the numerator scales like mp. The only terms that contribute through order g 2 m 3 D T are 



T 

J ■: 



(hs) 



3g+ig+gh 



g 2 T / A x (0,p 




+ 



-2g 2 - Ap 2 - Ad(n-Q) 2 + Am 2 D T Q ] A T (Q)A T (R) 
4r 2 - 2q 2 + 4p 2 l A T (Q)A X (R) 



-2dA T (Q) + 2^-(l-2q 2 A x (Q)] 



(46) 



In the terms that are already of order g 2 m 3 D T ', we can set R = —Q. In the terms of order 
g 2 rriDT s , we must expand the sum-integrand to second order in p. After averaging over 
angles of p, the linear terms in p vanish and quadratic terms of the form p % pP are replaced 
by p 2 8^ jd. We can set p 2 = —m 2 D , because any factor proportional to p 2 + m 2 D will cancel 
the denominator of the integral over p, leaving an integral with no scale. Our expression for 
the (hs) contribution reduces to 



J ■: 



(hs) 



3g+4g+gh 




+ N c g 2 m 2 D T 



Inserting the sum- integrals from Appendix [B] and the integrals from Appendix [C| this reduces 
to 

2 



(hs) 



■ig+Ag+gh 



37T 
11 

327T 



l m D T 3 



1 27 n 

— I h 27 

e 11 ' 



N c a s / n \ 2€ 
3?r UttT; 



2m D 



2c 



m 3 D T. 



(48) 



3. The (ss) contributions 



The (ss) contributions comes from the zero-frequency modes of the sum-integrals. The 
HTL correction functions Tp, T 000 , and T 0000 vanish when all the external frequencies are 

The 



m 



D 



l/(p 2 + m 2 D ) 



zero. The self-energy functions at zero-frequency are 11^(0, p) = and 11^(0, p) 
only scales in the integrals come from the longitudinal propagators Al(0, p) 
Therefore in dimensional regularization, at least one such propagator is required in order 
for the integral to be nonzero. The only terms in ( [34] ) that give nonzero contributions are 

N 

' ' ' ' ' A T (0,p)A r (0,q)A x (0,r) 



T-(ss) 

•^3g+4 S +5(/l 



4 



■9 2 T 2 



u 


-2q 2 - 4r 2 







+ 



p 2 + 4r 2 A T (0,p)A x (0,q)A x (0,r) 



(49) 
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After simplifying the integral by dropping terms that vanish in dimensional regularization, 
it reduces to 



T 

J ■ 



(ss) N c 2 2 r p 2 + Am 2 D 



g 2 T 2 f V + ^ m D (5Q) 

Jva v 2 (a 2 + m 2 „)(r 2 + m 2 n ) ' 



Wg+gh 4 y J pci p2( q 2 + m l)( r 2 + m l) 
Inserting the integrals from Appendix |C], this reduces to 



J ■: 



(88) 



3 



3g+4g+gh 



1 

- + 3 

e 



N c a s ( \i \ 4e 2 2 



37r \2mD 



m 2 D T 2 . (51) 



VI. THERMODYNAMIC POTENTIAL 

In this section, we calculate the thermodynamic potential Q(T, a s , mo, 5 = 1) explicitly, 
first to leading order in the 5 expansion and then to next-to-leading order. 

A. Leading order 

The complete expression for the leading order thermodynamic potential is the sum of 
the contributions from 1-loop diagrams and the leading term ( |19D in the vacuum energy 
counterterm. The contributions from the 1-loop diagrams, including all terms through order 
g 5 , is the sum of and (^0|): 

^i-loop = ^idcai < 1 - y«fl + S0m 3 D + — ( - + 2 log | - 7 + 2 7 + —J I , (52) 

where .Fideai is the free energy of an ideal gas of N 2 — 1 massless spin-one bosons, 

•Fideal = (N 2 - 1) (-^^ , (53) 

and mo and jl are dimensionless variables: 

mn = ^ , (54) 

A = ^- (55) 

Adding the counterterm (]T9|), we obtain the thermodynamic potential at leading order in 
the delta expansion: 

f 15 45 / u 7 vr 2 \ 1 

Olo = ^"ideai 1 - y«D + 30m^ + — f log I - - + T + T ^ ■ (56) 



The coefficient of m J, in fl5"6|) differs from the result calculated previously in Ref . [ 26 ] . In 



that paper, the constant under the logarithm of fi/2 was — |+7+log2 instead of — |+7+|7r 2 . 
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The reason for the difference is that the sum-integral T g was calculated in Ref. [26| using 
dimensional regularization to regularize the integral, but using the 3-dimensional expressions 
for the HTL propagators Ay and A^. At leading order, the difference can be absorbed 
into the definition of the scale /x. For calculations beyond leading order, it is essential for 
consistency to use the <i-dimensional expressions for these propagators.0 



B. Next-to-leading order 



The complete expression for the next-to-leading order correction to the thermodynamic 
potential is the sum of the contributions from the 2-loop diagrams, the HTL counterterms, 
and renormalization counterterms. The contributions from the 2-loop diagrams, including 
all terms through order g 5 , is the sum of (f4"5|), fl4"8|), and flSlj): 



n 



2— loop 



^ N c a s [ 15 165 

^ideal -= — \ - — + 45 m D — 

671 4 8 



I A 72 

- + 41og^ - — logm D + 1.969 



495 



I A 27 

- + 4 log £ - 2 log m D + — + 2 7 



D 



m 



D 



(57) 



The HTL counterterm contribution is the sum of (138|) and (1427): 



n 



HTL 




45 m 3 D -^Q + 21og^-7 + 2 7 



2tt 2 



Hi 



D 



(58) 



The ultraviolet divergences that remain after these 3 terms are added can be removed by 
renormalization of the vacuum energy density Sq and the HTL mass parameter mo- The 
renormalization contributions at first order in 5 are 



AQ = Ai£ + Aim|, 



d 



drrijj 



LO 



(59) 



where Ai£o and A\m 2 D are the terms of order 5 in the vacuum energy counterterm and 
the mass counterterm. The expression for A^o is given in (P7|). It cancels the poles in e 
proportional to in (^) and d58|). The remaining ultraviolet divergences are poles in e 
proportional to and m|j in d57|). If HTL perturbation theory (HTLpt) is renormalizable, 
both divergences must be removed by the same mass counterterm. This requires a remark- 
able coincidence between the coefficients of the two poles, and provides a nontrivial test of 
renormalizability. The value of the counterterm required is 



A l m 2 D 



11 N r a R 



4e 3tt 



m 



D 



(60) 



The complete contribution from the counterterms through first order in 5 is 



1 We thank E. Iancu and A. Rebhan for first bringing this problem to our attention. 
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(61) 



Adding the contributions from the two-loop diagrams in fl5"T|), the HTL counterterm in 
([58]) , and the renormalization counterterms in (|59|) and adding them to the leading order 
thermodynamic potential in ([56]), we obtain the complete expression for the thermodynamic 
potential at next-to- leading order in HTLpt: 



n 



NLO 



ideal ' 



1 - 15m D - 



N c a s 
3vr 



. p, 7 7T 2 

log 2-2 + 7 + y 



m 



D 



— + 45m D - — I log | - — logm D -2.001 J m£ 

495 / /i 5 \ 3 " ' 
+ — l0g 2 + 22 +7 ^ 



(62) 



C. Gap Equation 

The gap equation which determines mo is obtained by differentiating ( |6"2| ) with respect 
to mo and setting this derivative equal to zero yielding: 



m 



D 



i+ k*f- 5 + 7 + -g 



N c a s 
3vr 



11 / /t 36, \ 
— I log - - — logm D - 3.637 I m D 



5 

22 



(63) 



In Fig. |^, we have plotted the solution to this gap equation normalized to the leading-order 
perturbative result in (^) as a function of a s (2irT). The shaded band indicates the range 
resulting from varying the renormalization scale /i by a factor of two around /i = 2ttT. From 
this plot, we see that the gap equation solution matches nicely onto the perturbative result 
as a s — > 0. The solution decreases with a s (2nT) out to about a s ~ 0.06 and then begins to 
increase. It exceeds the perturbative result at around a s ~ 0.18, and then quickly diverges 
to +oo. 



VII. THERMODYNAMIC FUNCTIONS 



In this section, we compare the thermodynamic functions calculated at next-to-leading 
order in HTL perturbation theory (HTLpt) with those calculated using lattice gauge theory. 
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a s (27iT) 

FIG. 3. Solution to the gap equation ( |63[ ) as a function of a s (27rT). The shaded band corre- 
sponds to variation of the renormalization scale \x by a factor of two around \x = 2ttT. 



A. Pressure 



The final results for the LO and NLO HTLpt predictions for the free energy of pure- 
glue QCD are obtained by evaluating the thermodynamic potentials ( |56"D and (|62[) at the 
solution to the gap equation (|63D. Once the free energy JF(T) is given as a function of T, 
all other thermodynamic functions are determined. In particular, the pressure V and the 
energy density £ are 



V 



dT 



(64) 
(65) 



In Figure SI we have plotted the LO and NLO HTLpt predictions for the pressure of pure- 
glue QCD as a function of T/T c , where T c is the deconfinement transition temperature. To 
translate a s (27rT) into a value of T/T c , we use the two- loop running formula for pure-glue 
QCD with A^g = 0.65 T c . Thus a s (2nT) = 0.06 and 0.2 translate into T/T c = 415 and 
0.906, respectively. The LO and NLO HTLpt results are shown in Fig. |asa light-shaded 
band outlined by a dashed line and a dark-shaded band outlined by a solid line, respectively. 
The LO and NLO bands overlap all the way down to T = T c , and the bands are very narrow 
compared to the corresponding bands for the weak-coupling predictions in Fig. [I]. Thus the 
convergence of HTLpt seems to be dramatically improved over naive perturbation theory 
and the final result is extremely insensitive to the scale /x. 

In Fig. f|, w e have also included the 4-dimensional lattice gauge theory results of Boyd 
et al. || and the 3-dimensional lattice gauge theory results of Kajantie et al. [|5J. The LO 
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FIG. 4. The LO and NLO results for the pressure in HTLpt compared with 4-d lattice results 
(diamonds) and 3-d lattice results (dotted lines) for various values of an unknown coefficient in 
the 3-d effective Lagrangian. The LO HTLpt result is shown as a light-shaded band outlined by a 
dashed line. The NLO HTLpt result is shown as a dark-shaded band outlined by a solid line. The 
shaded bands correspond to variations of the renormalization scale \x by a factor of two around 
fi = 2vrT. 



and NLO HTLpt predictions differ significantly from the 4-d lattice results of Ref. H, even 
at the highest temperatures for which they are available. At T = 5 T c , the HTLpt prediction 
for the deviation of the pressure from that of the ideal gas is only 45% of the 4-d lattice 
result. In the high temperature limit, the HTLpt prediction approaches that of the ideal 
gas very slowly, in qualitative agreement with the results of the 3-d lattice calculations of 
Ref. M . However the quantitative agreement is not very good. The results of Ref. [7J depend 
on an unknown coefficient in the effective lagrangian for the dimensionally reduced theory. 
The 5 dotted lines in Fig. f| correspond to 5 possible values for that coefficient. We assume 
that the coefficient is such that the 3-d results match on reasonably well to the 4-d results, 
such as one of the middle 3 of the 5 dotted lines. In that case, the HTLpt prediction for the 
deviation from the ideal gas at T = 10 3 T C is only about 59% of the 3-d lattice result. We 
conclude that HTLpt at this order does not describe the pressure for pure- glue QCD. 



B. Trace Anomaly 



The combination 8 — 3V can be written as 

dT VT 4 



£-3?> = -T 5 4f^)- (CO 
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This combination is proportional to the trace of the energy-momentum tensor. In QCD with 
massless quarks, it is nonzero only because scale invariance is broken by renormalization 
effects. We will call it the trace anomaly density. It of course vanishes for an ideal gas of 
massless particles. However, it also vanishes for a gas of quasiparticles whose masses are 
linear in T and whose interactions are governed by a dimensionless coupling constant that 
is independent of T. 
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FIG. 5. The LO and NLO results for the trace anomaly in HTLpt. The LO HTLpt result 
is shown as a light-shaded band outlined by a dashed line. The NLO HTLpt result is shown as 
a dark-shaded band outlined by a solid line. The shaded bands correspond to variations of the 
renormalization scale /i by a factor of two around fj, = 2irT. 



In Fig. [|, we have plotted the LO and NLO HTLpt predictions for the trace anomaly 
density as a function of T/T c . At large T, the HTLpt prediction is very small and positive. 
As T decreases, the NLO prediction for £ — 3V increases to its maximum value around 10T C 
and then begins decreasing and quickly turns negative. The maximum value is less than 
about 0.2% of the energy density £ideai of the ideal gas. In contrast, the 4-d lattice result 
increases to a maximum of about 70% of £id ea i at a temperature that is very close to T c and 
then decreases rapidly to 0. 



VIII. CONCLUSIONS 



We have calculated the free energy of pure-glue QCD at high temperature to 2-loop or- 
der using HTL perturbation theory (HTLpt). The gauge invariance of the 2-loop expression 
was verified explicitly in generalized covariant gauge and generalized Coulomb gauge. The 
expression was reduced to a relatively compact form involving only scalar sum-integrals. 
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The numerical evaluation of the scalar sum-integrals would have been extremely difficult. 
We chose instead to approximate them by expanding in powers of m^/T, keeping all terms 
through 5 th order in g and m D /T. The ultraviolet divergences in the resulting expression 
for the thermodynamic potential can be removed by renormalization of the vacuum energy 
density and the HTL mass parameter m^. This provides a nontrivial test of the renormal- 
izability of HTL perturbation theory to this order. 

The two-loop order of HTLpt is the first order at which can be determined by a 
variational prescription. The condition that be a stationary point of the thermodynamic 
potential provides a "gap equation" for mo. The only ambiguity in the free energy then 
resides in the scale /z associated with renormalizations of the vacuum energy density and mr>. 
The predictions for the thermodynamic functions are extremely insensitive to the choice of 
fi. 

The quantitative predictions for the pressure in 2-loop HTLpt are disappointing. In the 
range 2T C < T < 20T C , the pressure is predicted to be nearly constant with a value of about 
95% of that of an ideal gas of gluons. The HTLpt prediction for the deviation from the 
ideal gas is about 45% of the result from 4-dimensional lattice gauge theory at T = 5T C , the 
highest temperature for which the lattice result is available. At very high temperature, the 
approach to the ideal gas limit is extremely slow, in qualitative agreement with the results 
of 3-d lattice gauge theory calculations. However, assuming that the 3-d results match on 
reasonably well to the 4-d results, the HTLpt prediction for the deviation from the ideal gas 
at T = 10 3 T c is only about 59% of the 3-d lattice result. 

A possible conclusion is that HTLpt at two-loop order is simply not a useful approxima- 
tion for thermal QCD. Another possibility is that the problem lies not with HTLpt but with 
our use of the mp/T expansion to approximate the scalar sum-integrals. The sum-integrals 
that were encountered at 4 th and 5 th order in m D /T were so difficult to evaluate that it 
seems hopeless to try to expand to higher order. However it is possible that the scalar 
sum-integrals could be evaluated numerically. Part of the difficulty is that it is necessary to 
isolate the infrared divergent and ultraviolet divergent terms analytically before evaluating 
the remaining terms numerically. Our m^/T expansions of the sum- integrals might be use- 
ful for generating the necessary subtractions that would allow the scalar sum-integrals to be 
evaluated numerically. 

Our calculations required the development of new methods for evaluating sum-integrals. 
The most difficult were two-loop sum-integrals that also involved an HTL angular average. 
These sum-integrals may be useful in other applications, such as solving the 2-loop 
derivable approximation for QCD. 
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APPENDIX A: HTL FEYNMAN RULES 



In this appendix, we present Feynman rules for HTL perturbation theory in pure-glue 
QCD. We give explicit expressions for the propagators and for the 3-particle and 4-particle 
vertices. The Feynman rules are given in Minkowski space to facilitate applications to real- 
time processes. A Minkowski momentum is denoted p = (po,p), and the inner product is 
p • q — po<?o ~ P • q- The vector that specifies the thermal rest frame is n = (1, 0). 



1. Gluon Self-energy 



The HTL gluon self-energy tensor for a gluon of momentum p is 

n^(p) = ml [T^(p, -p) - n"n"] . (A.l) 
The tensor T^ u (p, q), which is defined only for momenta that satisfy p + q = 0, is 

T^(p, -p) = (i/V— ) ■ (A.2) 

The angular brackets indicate averaging over the spatial directions of the light-like vector 
y = (l,y). The tensor T^ u is symmetric in p and v and satisfies the "Ward identity" 

p^(p,-p)=p-nn l/ . (A.3) 
The self-energy tensor IP" is therefore also symmetric in p and v and satisfies 

p^{p) = , (A.4) 
g^{p) = -m 2 D . (A.5) 

The gluon self-energy tensor can be expressed in terms of two scalar functions, the 
transverse and longitudinal self-energies n T and Ii L) defined by 

U L (p) = -H» , (A.7) 

where p is the unit vector in the direction of p. In terms of these functions, the self-energy 
tensor is 

TT{p) = -U T (p)T^ - \u L (p)I^ , (A.8) 



where the tensors T p and L p are 



V 


= <r- 


p»p v 


n^n u p 


TjXV 

L P 


n p n p 

n 2 p 







(A.9) 
(A.10) 
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The four-vector n£ is 



(A.ll) 



and satisfies p-n p = and n 2 = 1 — (n-p) /p . The equation ( |A.5| ) reduces to the identity 



(d-l)U T (p) + -U L (p) = m 2 D . 



nf 



(A.12) 



We can express both self-energy functions in terms of the function T 00 defined by ( |A.2 



n r (p) 



TO 



D 



T°°(p, -p) - 1 + n 2 



(d - l)n* 
= mfJl - T 00 (p, -p) 



(A.13) 
(A.14) 



In the tensor T M ^(p, — p) defined in ( |A.2| ), the angular brackets indicate the angular 
average over the unit vector y. In almost all previous work, the angular average in (|A.2|) 
has been taken in d = 3 dimensions. For consistency of higher order radiative corrections, it 
is essential to take the angular average in d = 3 — 2e dimensions and analytically continue 
to d = 3 only after all poles in e have been cancelled. Expressing the angular average as an 
integral over the cosine of an angle, the expression for the 00 component of the tensor is 



T>,-p) = ^/ 1 d C (i-c 2 r : 

Zi J — 1 i 



Po 



Po - P c 



(A.15) 



where the weight function w(e) is 



. . T(2 - 2e) n f 
w (e) = Arz f2 2e 



rm-eV ' r(|)r(i- e ) 



(A.16) 



The integral in (|A.15 ) must be defined so that it is analytic at p = oo. It then has a branch 
cut running from p = — |p| to p = +|p|. If we take the limit e — >• 0, it reduces to 



T>, -p) 



Po 
2|p| 



log 



Po + 


IpI 


Po - 


IpI 



(A.17) 



which is the expression that appears in the usual HTL self-energy functions. 

2. Gluon Propagator 

The Feynman rule for the gluon propagator is 

i5 ab A, u (p) , 



(A.18) 



where the gluon propagator tensor A^ v depends on the choice of gauge fixing. We consider 
two possibilities that introduce an arbitrary gauge parameter £ : general covariant gauge and 
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general Coulomb gauge. In both cases, the inverse propagator reduces in the limit £ — > oo 
to 



A^(pY v = -p 2 g^+p»p v - IP"(p) . 
This can also be written 



A T (p) p 



n 2 A L (p) p 



where A T and Ax are the transverse and longitudinal propagators: 

A r (p) ' 



Al(p) 



- n T (p) ' 
1 

-n£p 2 + ]J L (p) 



The inverse propagator for general £ is 

A-^r^A-^r-ipv 



A- 1 ( P r--(^-p-nn' 1 ) (f-v 



n n 



covariant 
Coulomb . 



The propagators obtained by inverting the tensors in ([A.24 ) and ( |A.23j ) are 



A""(p) = -A T (p)Tf + A L (p)nX P ~ £ 
= -A T (p)3f + A £ (p)nV-{ 



2\2 



(p 

(nip 1 



covariant , 
Coulomb . 



It is convenient to define the following combination of propagators: 

A x {p) = A L (p) + \a t {p) . 

Using ( A. 12 ), ( A.21| ), and ( |A.22| ), it can be expressed in the alternative form 



A*(p) 



m 



D 



dU T (p) A L {p)A T {p) 



(A.19) 
(A.20) 

(A.21) 
(A.22) 

(A.23) 
(A.24) 

(A.25) 
(A.26) 



(A.27) 



(A.28) 



which shows that it vanishes in the limit — > 0. In the covariant gauge, the propagator 
tensor can be written 



77 • T) 

A^(p) = [-A T (p)g^ + A x (p)n»n u ] - -^-A x (p) (p»n v + n»p v ) 



+ 



A T (p) + ^A x (p) - ^ 

pA pA 



pUpV 
pi 



(A.29) 



This decomposition of the propagator into three terms has proved to be particularly conve- 
nient for explicit calculations. For example, the first term satisfies the identity 

PnP X , n-p A x {p) 



[-A T (p) 9flu + A x (p)n„n„] A^(p) 



9p, 



p z 



+ 



n 2 p p 2 A L (p) 



P^p 



(A.30) 
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3. Three-gluon Vertex 



The three-gluon vertex for gluons with outgoing momenta p, q, and r, Lorentz indices 
H, u, and A, and color indices a, b, and c is 

KlciP, q, r) = -gf abc T^(p, q, r) , (A.31) 

where the three-gluon vertex tensor is 

T^\p, q, r) = gTip - qf + g v \q - rf + <? V (r - pf - m 2 D T^\p, q, r) . (A.32) 

The tensor T^ uX in the HTL correction term is defined only for p + q + r = 0: 

T^(p, q, r) = -/y"yY (-^ — ) ) . (A.33) 

This tensor is totally symmetric in its three indices and traceless in any pair of indices: 
g fMU T^ uX = 0. It is odd (even) under odd (even) permutations of the momenta p, q, and r. 
It satisfies the "Ward identity" 

q,T^\p, q, r) = T» x (p + q,r)- T x (p, r + q) . (A.34) 

The three-gluon vertex tensor therefore satisfies the Ward identity 

p^ x (p,q,r) = A-^-A-Vr*- (A.35) 



4. Four-gluon Vertex 

The four-gluon vertex for gluons with outgoing momenta p, q, r, and s, Lorentz indices 
/i, v, A, and a, and color indices a, b, c, and d is 

<r££T(p, q, r, s) = -ig 2 \j abx f xcd (g» x g^ - g»° g vX ) 

+2m^tr [T a (T b T c T d + T d T c T b )] T^ Xa (p, q, r, s) j 
+ 2 cyclic permutations , (A. 36) 

where the cyclic permutations are of (q,u,b), (r, A,c), and (s, <r, cf). The matrices T a 
are the fundamental representation of the SU (3) algebra with the standard normaliza- 
tion tr(T a T b ) = \5 ah . The tensor T^ vX(J in the HTL correction term is defined only for 
p + q + r + s = 0: 



T^(p,q,r,s) = (y^y x y c 



p-n 



K p-yq-y (q + r)-y 



(p + q)-n | (p + g + r)-n \ \ 



q-yr-y (r + s)-y r-y s-y (s + p)-y 
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This tensor is totally symmetric in its four indices and traceless in any pair of indices: 
g ^ v 7 '^ uXcr = 0. It is even under cyclic or anti-cyclic permutations of the momenta p, q, r, 
and s. It satisfies the "Ward identity" 

q^ vX °(p, q, r, s) = T vX °(p + q, r, s) - T» x "(p, r + q,s) (A.38) 

and the "Bianchi identity" 

T^ Xa (p, q, r, s) + T^ uXa (p, r, s, q) + T^ uXcT (p, s,q,r) = 0. (A.39) 

When its color indices are traced in pairs, the four-gluon vertex becomes particularly 
simple: 

5 ab 5 cd iT^(p, q, r, 8 ) = -ig 2 N c (N 2 c - l)T^(p, q, r, a) , (A.40) 

where the color-traced four-gluon vertex tensor is 

T^(p, q, r, s) = 2^V a - 9*ST - gT<r X - m 2 D T^(p, s, q, r) . (A.41) 

Note the ordering of the momenta in the arguments of the tensor T^ vXa , which comes from 
the use of the Bianchi identity ( |A.39| ). The tensor ( [A.41]) is symmetric under the interchange 
of \x and u, under the interchange of A and a, and under the interchange of (/i, v) and (A, a). 
It is also symmetric under the interchange of p and q, under the interchange of r and s, and 
under the interchange of (p, q) and (r, s). It satisfies the Ward identity 

p^ x °(p, q, r, s) = T uX °(q, r + p,s)- T uX °(q, r jS + p). (A.42) 



5. Ghost propagator and vertex 

The ghost propagator and the ghost-gluon vertex depend on the gauge. The Feynman 
rule for the ghost propagator is 

^r5 ab covariant , (A.43) 

p l 

-^5 ab Coulomb . (A.44) 

n l v p l 

The Feynman rule for the vertex in which a gluon with indices fi and a interacts with an 
outgoing ghost with outgoing momentum r and color index c is 

_ g fabc r fi covariant , (A.45) 

_ g fabc ^_ r . nn ^ Coulomb . ( A.46) 

Every closed ghost loop requires a multiplicative factor of —1. 
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6. HTL Counterterm 



The Feynman rule for the insertion of an HTL counterterm into a gluon propagator is 

-i5 ab W v (p), (A.47) 
where IP^p) is the HTL gluon self-energy tensor given in ( |A.8|) . 

7. Imaginary-time formalism 

In the imaginary-time formalism, Minkoswski energies have discrete imaginary values 
Po = i(2TtnT) and integrals over Minkowski space are replaced by sum-integrals over Eu- 
clidean vectors (27rnT, p). We will use the notation P = (Pq,p) for Euclidean momenta. 
The magnitude of the spatial momentum will be denoted p = |p|, and should not be confused 
with a Minkowski vector. The inner product of two Euclidean vectors is P ■ Q = PqQq + p ■ q. 
The vector that specifies the thermal rest frame remains n = (1, 0). 

The Feynman rules for Minkowski space given above can be easily adapted to Euclidean 
space. The Euclidean tensor in a given Feynman rule is obtained from the corresponding 
Minkowski tensor with raised indices by replacing each Minkowski energy p by iPo, where 
Pq is the corresponding Euclidean energy, and multipying by —i for every index. This 
prescription transforms p = (po, p) into P = (Pq, p), g^ v into —5^, and p-q into —P-Q. The 
effect on the HTL tensors defined in (|A.2|) , (|A.33|) , and (|A.37|) is equivalent to substituting 



p-n — ► —P-N where N = (— i, 0), p-y — > — P-Y where Y = (— i, y), and — > Y 1 *. For 
example, the Euclidean tensor corresponding to ( |A.2j ) is 



/ P-N\ 

r^(P, -P) = \ rtx Y u —— \ . (A.48) 

The average is taken over the directions of the unit vector y. 

Alternatively, one can calculate a diagram by using the Feynman rules for Minkowski 
momenta, reducing the expressions for diagrams to scalars, and then make the appropriate 
substitutions, such as p 2 — > —P 2 , p ■ q — > — P ■ Q, and n ■ p — > in ■ P. For example, the 
propagator functions ( A..21 ) and ( A.22|) become 



A - (/> > = wrkm ' (A49) 
A '( p > = wrkm ■ (A ' 50) 



The expressions for the HTL self-energy functions Ht(P) and Hl(P) are given by (|A.13|) 
and (|A.14|) with n 2 replaced by rip = p 2 /P 2 and T 00 (p, —p) replaced by 



rp _!^£* & _^_£_. (A . 51) 

Note that this function differs by a sign from the 00 component T 00 (P, — P) of the Euclidean 
tensor corresponding to (|A.2p : 
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= -T P . (A.52) 

po—>iPo 



A more convenient form for calculating sum-integrals that involve the function Tp is 



2 



Tp = ( JfT^jt • (A ' 53) 



where the angular brackets represent an average over c defined by 



(/(c)) c = «,( e ) dc (1 - c 2 )~7(c) (A.54) 
Jo 



and w(e) is given in ( [A.16Q . 



APPENDIX B: SUM-INTEGRALS 

In the imaginary-time formalism for thermal field theory, a boson has Euclidean 4- 
momentum P = (Pq, p), with P 2 = P 2 + p 2 . The Euclidean energy Pq has discrete values: 
Po = 2irnT, where n is an integer. Loop diagrams involve sums over Pq and integrals over p. 
With dimensional regularization, the integral is generalized to d = 3 — 2e spatial dimensions. 
We define the dimensionally regularized sum-integral by 

where d = 3 — 2e is the dimension of space and \x is an arbitrary momentum scale. The factor 
(e 7 /47r) e is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, \x coincides with the renormalization scale of the MS renormalization scheme. 



1. Simple one-loop sum-integrals 



The simple one-loop sum-integrals required in our calculations are 



_!Lt 4 

45 ' 

WJM 2e ± 

\AttTJ 12 



logP 2 
1 

P (P 2 ) 2 

ip T P 2 ~ (4vr) 2 V4vrT. 
r 1 1 /^x* 

/p (P 2 ) 2 (4tt) 2 Uttt; 



i / ^ 2e 



1 / 7T 2 \ 

- + 2 7 + 2 + I 4 + 4 7 + — - 471 j e 



i + 2 T + f ^ - 4 7 i j e 



(B.2) 
(B.3) 
(B.4) 
(B.5) 
(B.6) 
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p 



1 / fj, 



p (P 2 ) 3 

4 



2c 



P 



(4tt) 2 \4irTj 4 
1 / a \ 2e 5 
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p (P 2 ) 4 
1 

2Y3 



p (P 



(4tt) 2 V4ttT 
2C(3) 1 
(4vr) 4 T 2 ' 



1 16 
- + 2 7 -- 
e 15. 



(B.7) 
(B.8) 
(B.9) 



The calculation of these sum-integrals is standard. The errors are all one order higher in e 
than the smallest term shown. The number 71 is the first Stieltjes gamma constant defined 
by the equation 



C(l + *) = - + 7 - liz + 0{z 2 ) 



(B.10) 



2. One-loop HTL sum-integrals 

The one-loop sum-integrals involving the HTL function Tp defined in ( |A.51| ) are 



pP^ Tp 



Tp 



T Gttt) ( 24) 
1 / y, ^ 2e 



1 c'(-i) 

e C-l 



p ;r 
1 

p p 2 P 2 
1 



p (P 2 ) 2 
1 

p p 4 



(4tt; 


2 


UvrT 


1 




/ /i 


(4vr; 




UvrT 


1 




f " 


(4< 




UttT 


1 




f " 


(4vr; 




UttT 



-1 



+ 2 7 + 2 log 2 



2c 



2c 



2c 



I \ 7T 2 

21og2(- + 2 7 J +21og 2 2 + y 



- + 2 7 + 1 

e 



l + 21og2) (- + 2 7 



^ + ylog2 + 21og 2 2 



(B.ll) 
(B.12) 
(B.13) 
(B.14) 

(B.15) 



The errors are all of order e. 

It is straightforward to calculate the sum-integrals (|B.11|) -( |B.15|) using the representa- 
tion ( |A.53j ) of the function Tp. For example, the sum-integral ( |B.11| ) can be written 



tp P 2 



T P 



p2 



pP 2 +p 2 \P 2 +p 2 c 2 , 



(B.16) 



of P 2 in the numerator to cancel denominators, this becomes 



where the angular brackets denote an average over c as defined in ( A.54|) . Using the factor 

(B.17) 



p T* Tp 



1-c 2 Tp\P 2 P 2 +p 2 c 2 
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After rescaling the momentum by p — > p/c, the second sum- integral on the right side 
becomes the same as the first sum-integral, and the expression reduces to 



p P z 



-,-1+26 




(B.18) 



Evaluating the average over c, using the expression ( |B.3|) for the sum-integral, and expanding 
in powers of e, we obtain the result ( |B.11|) . Following the same strategy, all the sum-integrals 
( |B.11| )-( |B.15| ) can be reduced to linear combinations of the simple sum-integrals ( |B.3| ) and 



( p.5|) with coefficients that are averages over c. The only difficult integral is the double 
average over c that arises from (B.15): 



' 3+2e 3+2e ' 
Ci — C 2 



l + 21og2 



r 2 



' ci,c 2 



9 9 5 3 6 J 



(B.19) 



3. Simple two-loop sum-integrals 



The simple two-loop sum-integrals that are needed are 
1 



pq P 2 Q 2 R 2 
1 







r- 



pq P 2 Q 2 r 2 
9 2 



m i 

(4vr) 2 UvrT/ 12 



i + 10 

e 



121og2 + 4 



C'(-i) 
C(-i)J 



T' 1 



pq P 2 Q 2 r 4 
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pq P 2 Q 2 r 2 R 2 

PQ 



MM 



k 



t 



(An) 2 \4ttT ) 6 
T 2 / fi \ 4e 1 
(4tt) 2 V4^T/ 9 

^2 / „ \4e 



e 3 + 7+ C-l 



- + 7.521 

e 



pq P 2 Q 2 r 4 (An) 



{atit) ( 8) 



4C 



12 8 

— I h 4 log 2 + -7 H 

e ^ 9 ^ logz ^ 3 7 3 C -l) 



(B.20) 
(B.21) 
(B.22) 
(B.23) 
(B.24) 



where R = —(P + Q) and r = |p + q|. The errors are all of order e. 

To motivate the integration formula we will use to evaluate the two-loop sum-integrals, 
we first present the analogous integration formula for one-loop sum-integrals. In a one-loop 
sum-integral, the sum over Pq can be replaced by a contour integral in po = —iPq: 



4 F(P) = lim / [F(-zp ,p) -F(0,p)]e^°n(p 

Tp n->o+ J 2m J-p 



(B.25) 



where n(po) = l/(e /3po — 1) is the Bose-Einstein thermal distribution and the contour runs 
from — oo to +oo above the real axis and from +oo to — oo below the real axis. This formula 
can be expressed in a more convenient form by collapsing the contour onto the real axis and 
separating out those terms with the exponential convergence factor n(|p |)- The remaining 
terms run along contours from — oo ± ie to and have the convergence factor e vpo . This 
allows the contours to be deformed so that they run from to ±ioo along the imaginary 
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Po axis, which corresponds to real values of Pq = —ipo- Assuming that F(— ipo,p) is a real 
function of p , i.e. that it satisfies F(— ipo,p) = F(—ip ,p)*, the resulting formula for the 
sum-integral is 

f p F{P) = j p F{P) + Je(p )n(\ Po \)2hnF(-i Po + e,p) , (B.26) 

where e(po) is the sign of po- The first integral on the right side is over the (d+l)-dimensional 
Euclidean vector P = (Po, p) and the second is over the (d + l)-dimensional Minkowskian 
vector p — (p , p). 

The two-loop sum-integrals can be evaluated by using a generalization of the one-loop 
formula ( gg§) : 

4 F(P)G(Q)H(R) = \f F(P)G(Q)H(R) 

JPQ 3 JPQ 

+ [ e(p )n(\p \)2ImF(-ip + e,p)Re [ G(Q)H(R) 

Jp J Q P =-ipo+£ 

+ e(po)n(\p Q \)2ImF(-ip + e,p) / e(q )n(\q \) 2lmG(-iq + e, q) 

xReif( J R) + (cyclic permutations of F, G, H) . (B.27) 

Bo=i(Po+9o)+e 

The sum over cyclic permutations multiplies the first term on the right side by 3, so there 
are a total of 7 terms. This formula can be derived in 3 steps. First, express the sum over 
Po as the sum of two contour integrals over po, one that encloses the real axis Impo = and 
another that encloses the line Impo — — Im<7o- Second, express the sum over go as a contour 
integral that encloses the real-go axis. Third, symmetrize the resulting expression under the 
6 permutations of F, G, and H. The resulting terms can be combined into the expression 
( |B.27p . The integrals of the imaginary parts that enter into our calculation can be reduced 
to 



/ e(p )n(M) 2Im^ f(-ip + e, p) = / ^ \ £ f(±ip + s,p), (B.28) 

e(po)n(\p \) 2lmT P 

Pa=-ipo+e 

= ~J pn(p)^( c ' 3+2€ f( ± W + s,p/c)) c . (B.29) 



P =-ip +£ JP P * ± 

f(-ip + s,p) 

Po=-ipo+£ 



The latter equation is obtained by inserting the expression ( |A.53| ) for 7p, using ( p.28| ), and 
then making the change of variable p — > p/c to put the thermal integral into a standard 
form. 

As a simple illustration, we apply the formula ( B.27| ) to the sum-integral ( B.21 ). The 
nonvanishing terms are 



'pq ^gv = 2 27r6{pl ~ p2) J Q oh 

+ f n(\p \) 2<K8(pl - p 2 ) f n(\q \) 2n5(q 2 - q 2 )\ . (B.30) 

Jp J a T 
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The delta functions can be used to evaluate the integrals over po and go- The integral over 
Q is given in ( |C.98|) up to corrections of order e. This reduces the sum-integral to 



pq P 2 Q 2 r 2 (4vr) 



4-21og2 



2c 



U (P) p -2e + 
P /' ■' pq 



n(p)n(q) 1 



pq 



(B.31) 



The momentum integrals are evaluated in ( |C.5| ) and flC.6j ). Keeping all terms that contribute 
through order e°, we get the result ( |B.21| ). The sum-integral ( |B.22| ) can be evaluated in the 
same way: 



Q 



pq P 2 Q 2 r i (4vr) 



-21og2 



2c 



<P) p -2c + 
p P -/pq 



n(p)n(q) q 2 
pq r 4 



[B.32) 



The sum-integral (P-24Q can be reduced to a linear combination of ( |B.21|) and ( |B.22| ) by 



expressing the numerator in the form P-Q = PqQo + {r 2 — p 2 — q 2 )/2 and noting that the 
PoQo term vanishes upon summing over Po or Qq. 

The sum-integral (|B.23| ) is a little more difficult. After applying the formula QB.27] ) and 
using the delta functions to integrate over p , q , and r , it can be reduced to 



n(p) 



pq P 2 Q 2 r 2 R 2 J P p Jq Q 2 R 2 \r 
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pq 
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+ -z + 



q 2 q 2 I A(p,q,r) 



(B.33) 



where A(p,q,r) is the triangle function that is negative when p, q, and r are the lengths of 
3 sides of a triangle: 



A(p, q,r) = p +q +r —2(pq +qr +rp). 



(B.34) 



After using (|C.104|) -( |C.106|) to integrate over Q, the first term on the right side of (B.33| ) is 
evaluated using ( |C.5|) . The 2-loop thermal integrals on the right side of ( p.33|) are given in 
( |C.8|) - (|C.11| ). Adding together all the terms, we get the final result (|B.23|) . 



4. Two-loop HTL sum-integrals 



The two-loop sum-integrals involving the HTL function Tp defined in ( |A.51|) are 
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(B.35) 



^B.36) 
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P-Q 



] p Q P 2 Q 2 r^ R (4vr) 2 UtttJ 



T 2 ( fi \ 4e 

~96 
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To calculate the sum- integral (|B.35|) , we begin by using the representation ( A.53|) of the 
function T R : 
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The first sum- integral on the right side is given by (|B.21| ). To evaluate the second sum- 
integral, we apply the sum- integral formula (|B.27|) : 
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where r c = |p + q/c|. In the terms on the right side with a single thermal integral, the 
appropriate averages over c of the integrals over Q are given in ( CTTU5 ) and ( |C.102| ) . 
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The subsequent integral over p is a special case of (|C.5[) : 
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where (a)& = T(a + b)/T(a) is Pochhammer's symbol. Combining this with ( [B.40| ), we obtain 
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For the two terms in ( |B.39|) with a double thermal integral, the averages weighted by c 2 are 
given in ( |C.12| ) and (|C.15| ). Adding them to (|B.42|) , the final result is 



k 



pq P 2 Q 2 \B% + r 2 c 2 I (4vr) 2 \AttTJ 48 



^ - 12 log 2 + 4|p^ ~ + 18.66 



(B.43) 



35 



Inserting this into ( |B.38| ), we obtain the final result ( |B.35| ). 

The sum- integral (|B.36| ) is evaluated in a similar way to ( |B.35|) . Using the representation 
for T R , we get 
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The first sum-integral on the right hand side is given by ( B.22 ). To evaluate the second 
sum-integral, we apply the sum- integral formula ( B.27Q : 
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In the terms on the right side with a single thermal integral, the weighted averages over c 
of the integrals over Q are given in ( |C.112j ), ( |C.113| ), and (|C.108| ): 
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After using ( [B.41 ) to evaluate the thermal integral, we obtain 
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For the two terms in ( |B.45|) with a double thermal integral, the averages weighted by c 2 are 
given in ( |Cl4T ), ([0171) , and ( |Cl8| ). Adding them to (|B~47|) , the final result is 
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Inserting this into (|B.44|) , we obtain the final result (|B.36| ). 

To evaluate ( |B.37| ), we use the expression ( |A.53| ) for Tr and the identity P-Q = (i? 2 
P 2 — Q 2 )/2 to write it in the form 
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The sum-integrals in the first 3 terms on the right side of ( p.49|) are given in ( |B.3|) , flB.12|) , 
( |B.21 ), and QB.24Q . The last sum- integral before the average weighted by c is given in flB.38 ). 
The average weighted by c 2 is given in ( |B.43| ). The average weighted by c 4 can be computed 
in the same way. In the integrand of the single thermal integral, the weighted averages over 
c of the integrals over Q are given in ( |C.111| ) and ( |(J.103| ) : 
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After using ( |B.41| ) to evaluate the thermal integral, we obtain 
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For the two terms with a double thermal integral, the averages weighted by c 4 are given in 
(pT3|) and ( |Ol6|) . Adding them to flB31|) , we obtain 
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Inserting this into (|B.49|) along with (|B.43|), we get the final result (|B.37|). 



APPENDIX C: INTEGRALS 

Dimensional regularization can be used to regularize both the ultraviolet divergences 
and infrared divergences in 3-dimensional integrals over momenta. The spatial dimension is 
generalized to d = 3 — 2e dimensions. Integrals are evaluated at a value of d for which they 
converge and then analytically continued to d — 3. We use the integration measure 

_ /e> 2 Y f d^p 



1. 3-dimensional integrals 



We require several integrals that do not involve the Bose-Einstein distribution function. 
The momentum scale in these integrals is set by the mass parameter mo. The one-loop 
integrals are 
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We also require a two-loop integral: 
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The errors in (|C.2|) - (|C.4|) are all one order higher in e than the smallest term shown. 

2. Thermal integrals 



The thermal integrals involve the Bose-Einstein distribution n(p) = \j{ef p 
one-loop integrals can all be obtained from the general formula 

f n{pl = C(2 + 2 a -2e)r(2 + 2a-26)r(|) , y 
hp 4vr 2 r(| - e) V M / 

The simple two-loop thermal integrals that we need are 
r n(p)n(q) 1 T 2 / y, \ 4e / 1 
ipq pq r 2 ~ (4vr) 2 V4^T/ V 4 
n(p)n(q) p 2 T 2 f \x \ 4f 



2a-2e 



1 14 , C'(-l) 



pq 



pq 



" 4 (4tt) 2 UttT/ 



9 3 ' 3 C(-l) 



4.855 e 



1). The 
(C.5) 

(C.6) 
(C.7) 



We also need some more complicated 2-loop thermal integrals that involve the triangle 
function defined in (|B.34|) : 
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The most difficult thermal integrals to evaluate involve both the triangle function and the 
HTL average defined in (|A.54 ). There are 2 sets of these integrals. The first set is 
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The second set of these integrals involve the variable r c = |p + q/c|: 
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The errors in (|C.6|) - (|C.1§ ) are all one order higher in e than the smallest term shown. The 
numerical constant in ( |U.8| ) can be expressed analytically in terms of the transcendental 
numbers appearing in (|C.9| ) and ( |C.10|) . We do not know how to calculate the numerical 
constants in QC77D , ( PTT| ), (P^2])-(P^) analytically. 

The simplest way to evaluate integrals like ( C.6 ) and flCTf ) whose integrands factor into 
separate functions of p, q, and r is to Fourier transform to coordinate space where they 
reduce to an integral over a single coordinate R: 
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The Fourier transform is 
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and the dimensionally regularized coordinate integral is 
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The Fourier transforms we need are 
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If a is an even integer, the Fourier transform ( |C.23| ) is particularly simple in the limit d — > 3: 
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where x = irRT. We can use these simple expressions only if the integral over the coordinate 
R in ( p.!9|) converges for d = 3. Otherwise, we must first make subtractions on the integrand 
to make the integral convergent. 

The integral (|C.7|) can be evaluated directly by applying the Fourier transform formula 
7T9|) in the limit e — > 0. The integral (p.6|) however requires subtractions. It can be 
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In the second term on the right side, the integral over q is proportional to p _1_2e , so the 
integral over p can be evaluated using ( |C5| ). This first term on the right side is convergent 
for d = 3 so it can be evaluated easily using the Fourier transform formula ( |C.19|) . The 
integral over R reduces to a sum of integrals of the form Jq°° dx x m coth n x. Although the 
sum of the integrals converges, each of the individual integrals diverges either as x — > or 
as x — > oo. A convenient way to evaluate these integrals is to use the strategy in Appendix 
C of Ref. M . The integrals are regularized by using the substitution 



r 00 T(l + 8) 

/ dxx m coth n x ' 
Jo 2 s 



dx x m coth n x — ► 1 T ; / dx x m+s coth n x . (C.27) 



The divergences appear as poles in 5 that cancel upon adding a convergent combination of 
these integrals. 

The integrals ( |C8[ )-( P?TD| ) can be evaluated by first averaging over angles. The triangle 
function can be expressed as 

A(p,q,r) = -4p 2 q 2 (l-cos 2 6) , (C.28) 
where 9 is the angle between p and q. For example, the angle average for ( |(J.8| ) is 
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After integrating over x and inserting the result into ( |C.8| ), the integral reduces to 
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The integrals over p and q factor into separate integrals that can be evaluated using 
After averaging over angles, the integrals (|C.9|) and flC.10|) reduce to 
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The integral (|C.11|) can be evaluated by using the remarkable identity 
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The identity can be proved by expressing the angular averages in terms of integrals over the 
cosine of the angle between p and q as in ( |(J.29j ), and then integrating by parts. Inserting 
the identity (|C.33| ) into ( |C.11| ), the integral reduces to 
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The integral in the first term on the right is given in flCTf ), while the second term can be 
evaluated using ( |C.5|) . 

To evaluate the weighted averages over c of the thermal integrals in (|C.12j )- (|C.14D , we 
first isolate the divergent parts, which come from the region p — q — ► 0. We write the product 
of thermal functions in the form 



n(p)n(q) = n(p)n(q) — 



s 2 n 2 (s) 



pq 



+ 



s 2 n 2 (s) 
pq 



(C.35) 



where s = (p + q)/2. In the difference term, the HTL average over c and the angular average 
over x = p • q can be calculated in 3 dimensions: 
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The remaining 2-dimensional integral over p and q can then be evaluated numerically: 
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The integrals involving the n 2 (s) term in (|C.35|) are divergent, so the HTL average over c 
and the angular average over x = p • q must be calculated in 3 — 2e dimensions. The first 
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step in the calculation of the n 2 (s) term is to change variables from p and q to s = (p + q)/2, 
(3 = 4pq/ (p + q) 2 , and x = p • q: 
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x, 2 f d(3 r 2e (l - P)- 1/2 (f(s + , s-,r) + f(s-,s+, r)\ , (C.42) 



where s± = s[l ± y/1 — (3\ and r = s[4 — 2/3(1 — x)] 1//2 . The 2 terms inside the average over x 
come from the regions p > q and p < q, respectively. The integral over s is easily evaluated: 
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It remains only to evaluate the averages over c and x and the integral over j3. 

The first step in the calculation of the n 2 {s) term of ( C.12 ) is to decompose the integrand 
into 2 terms: 
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The weighted averages over c give hypergeometric functions: 
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In the +q case of (|C.45| ), the ie prescription is unnecessary. The argument of the hyper- 
geometric function can be written 1 — (3y, where y = (1 — x)/2. After using a transformation 
formula to change the argument to f3y, we can evaluate the angular average over x to obtain 
hypergeometric functions with argument (3. For example, the average over x of (|C.45|) is 
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where (a)b is Pochhammer's symbol which is defined in ( |U.127| ). Integrating over (3, we 
obtain hypergeometric functions with argument 1: 
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The integral weighted by c 4 can be evaluated in a similar way. Expanding in powers of e, 
we obtain 
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In the — q case of ( |C.45| ), the argument of the hypergeometric functions can be written 
(1— (3y)/(l— /3±ze), where y = (l—x)/2 and the prescriptions +ie and — is correspond to the 
regions p > q and p < q, respectively. These regions correspond to the two terms inside the 
average over x in (|C.42|) . In order to obtain an analytic result in terms of hypergeometric 
functions, it is necessary to integrate over (3 before averaging over x. The integrals over 
(3 can be evaluated by first using a transformation formula to change the argument of the 
hypergeometric function to —/3(l — y)/(l — (3) and then using the integration formula ( |(J.134| ) 
to obtain hypergeometric functions with arguments y or 1 — y: 
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After averaging over x, we obtain hypergeometric functions with argument 1: 
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The integral weighted by c 4 can be evaluated in a similar way. Expanding in powers of e 
and then taking the real parts, we obtain 
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Inserting the sum of the integrals (p.49|) and ( |C.53|) into the thermal integral ( |C.42| ) and 
similarly for the integrals weighted by c 4 , we obtain 
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Adding these integrals to the subtracted integrals in ( |C.39|) - (|C.40| ), we obtain the final 
results in ( |CH3 )-(|CHp. 



To evaluate the subtraction in the integral ( |C.41|) , we use the identity q 2 = (r 2 + q 2 — 
p 2 — 2p • q)/2. The integral with q 2 — p 2 in the numerator is purely imaginary. Thus the 
real part of the integral can be expressed as 
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It remains only to evaluate the integral with p • q in the numerator. We begin by using the 
identity 
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In the first term on the right side, the average over c is a simple multiplicative factor: 
(c 2 ) c = 1/(3 — 2e). The average over x gives hypergeometric functions of argument (3: 
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The integral over (3 gives hypergeometric functions of argument 1: 
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Expanding in powers of e, we obtain 
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In the second term of ( U.58 ), the average over c is given by ( |C.46| ). In the +q term, the 
average over x = p • q is 
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Integrating over j3, we obtain hypergeometric functions of argument 1: 
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Expanding in powers of e, we obtain 
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In the — q term in the integral of the second term of ( |C.58|) , we integrate over f3 before 
averaging over x. The integral over f3 can be expressed in terms of hypergeometric functions 
of type 2-F1: 
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The phase in the last term is e _Mre for the f(s + , s_, r) term of ( |C.42 ), which comes from the 
p > q region of the integral, and e t7Te for the /(s_, s+, r) term, which comes from the p < q 
region. The average over x = p • q can be expressed in terms of hypergeometric functions 
of type 3F2 evaluated at 1: 
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The expansion of the real part of the integral in powers of e is 
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Inserting ( |C.61|) , (|C.64|) , and ( |C.67| ) into the thermal integral of ( |C.58|) , we obtain 
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Inserting this along with ( |(J.55| ) into ( |U.57| ), we obtain 
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Adding this integral to the subtracted integral in ( C.41 ), we obtain the final result in (|C.14 ). 

To evaluate the weighted averages over c of the thermal integrals in ( |CJ . 1 5| )— ( |CJ ■ 1 8| ) , we 
first isolate the divergent parts, which arise from the region q — > 0. For the integrals ( |C.15| ) 
and (|C.16j ), a single subtraction of the thermal distribution n(q) suffices to remove the 
divergences: 
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For the integral (|C.17|) , a second subtraction is also needed to remove the divergences 
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In the last integral ( |C.18 ), it is convenient to first use the identity r 2 = p 2 + 2p ■ q/c+ q 2 /c 2 
to expand it into 3 integrals, two of which are ( |C15| ) and ( |C17| ). In the third integral, the 
subtraction (|C.71| ) is needed to remove the divergences. For the convergent terms, the HTL 
average over c and the angular average over x = p ■ q can be calculated in 3 dimensions: 
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The remaining 2- dimensional integral over p and q can then be evaluated numerically: 
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The integrals involving the terms subtracted from n(q) in ( |U.70 ) and ( |C.71| ) are divergent, 
so the HTL average over c and the angular average over x = p ■ q must be calculated in 
3 — 2e dimensions. The first step in the calculation of the subtracted terms is to replace the 
average over c of the integral over q by an average over c and x: 
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The integral over p can now be evaluated easily using either (|B.41| ) or 
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It remains only to calculate the averages over c and x. The averages over x give 2 -Pi hyper- 
geometric functions with argument [(1 =F c)/2 — ie]^ 1 : 
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Using a transformation formula, the arguments can be changed to (1 =F c)/2 — is. If the 
expressions ( p.81| ) and ( |CJ.82| ) are averaged over c with a weight that is an even function of 
c, the + and — terms combine to give 3 F 2 hypergeometric functions with argument 1. For 
example, 
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Upon expanding the hypergeometric functions in powers of e and taking the real parts, we 
obtain 
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If the expressions ( p.82|) and (|C.82j ) are averaged over c with a weight that is an odd function 
of c, they reduce to integrals of 2 F\ hypergeometric functions with argument y. For example, 
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The expansions of the integrals of the hypergeometric functions in powers of e are given in 
( |C.147|) -( |C.148|) . The resulting expansions for the real parts of the averages over c and x are 
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Multiplying each of these expansions by the appropriate factors from the integral over q in 
( p.79| ) and the integral over p in ( |C.80| ) or ( |B.41| ), we obtain 
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Adding these integrals to the subtracted integrals in ( |C.75|) - (|C.77| ), we obtain the final 
results in dC15D - (|CT?D . Combining (|CT78|) with flC^l and ( P^96|) , we obtain 
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The final integral (|C.18|) is obtained from ( |C.15|) , ( |C.17| ), and ( |C.97|) by using the identity 
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3. 4-dimensional integrals 



In the sum-integral formula (|B.27|) , the second term on the right side involves an integral 
over 4-dimensional Euclidean momenta. The integrands are functions of the integration 
variable Q and R — —(P+Q). The simplest integrals to evaluate are those whose integrands 
are independent of Pq : 
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Another simple integral that is needed depends only on P 2 = P 2 + p 2 : 
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where (a)b is Pochhammer's symbol which is defined in ( |U.127| ). We need the following 
weighted averages over c of this function evaluated at P = (—ip, p/c): 
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The remaining integrals are functions of Pq that must be analytically continued to the 
point P = —ip + e. Several of these integrals are straightforward to evaluate: 
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We also need a weighted average over c of the integral in (|C.104|) evaluated at P 
The integral itself is 
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The weighted average is 
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The most difficult 4- dimensional integrals to evaluate involve an HTL average of an 
integral with denominator Rq + r 2 c 2 : 
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The analytic continuation to Po = ~~ W + £ is implied in these integrals and in all the 
4-dimensional integrals in the remainder of this subsection. 

We proceed to describe the evaluation of the integrals ( |C.109| ) and flClllj ). The integral 
over Q Q can be evaluated by introducing a Feynman parameter to combine Q 2 and _Rq + r 2 c 2 
into a single denominator: 
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where we have carried out the analytic continuation to P — ~W + £ - Integrating over r 
and then over the Feynman parameter, we get a 2 -Pi hypergeometric function with argument 
1-c 2 : 
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The subsequent weighted averages over c give 3F2 hypergeometric functions with argument 
1: 



1 
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(C.116) 
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(C.117) 



After expanding in powers of e, the real part is ( U.lll] ). 

The integral ( C.112j ) has a factor of 1/r 2 in the integrand. After using (|C.114|) , it is 
convenient to use a second Feynman parameter to combine (1 — x + xc 2 )r 2 with the other 
denominator before integrating over r: 



Q 2 r 2 (Rl + r 2 c 2 ) 
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3 f 1 f 1 
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8 Jo Jo 



1 — x + xc )r + 2y(l — x)r-p + y(l — x) p — is 
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(C.118) 



After integrating over r and then y, we obtain 2 F\ hypergeometric functions with arguments 
x{l — c 2 ). The integral over x gives a 2-F1 hypergeometric function with argument 1 — c 2 : 
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After averaging over c, we get a hypergeometric function with argument 1: 
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(C.120) 



After expanding in powers of e, the real part is ( U.112j ). 

To evaluate the integral ( C.113|) , it is convenient to first express it as the sum of 3 
integrals by expanding the factor of q 2 in the numerator as q 2 = p 2 + 2p ■ r + r 2 : 



Q 2 r 2 (Rl + r 2 c 2 ) Jq \r 



P +2^ + 1 



Q 2 (R 2 + r 2 c 2 ) 



(C.121) 



To evaluate the integral with p ■ r in the numerator, we first combine the denominators 
using Feynman parameters as in (|C.11S| ). After integrating over r and then y, we obtain 
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2-F1 hypergeometric functions with arguments x(l — c 2 ). The integral over x gives 2^1 
hypergeometric functions with arguments 1 — c 2 : 
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After averaging over c, we get a hypergeometric function with argument 1: 
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After expanding in powers of e, the real part is 
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Combining this with ( |C.109|) and flC.lll ), we obtain the integral ( C.113|) . 



4. Hypergeometric functions 



The generalized hypergeometric function of type p F q is an analytic function of one vari- 
able with p + q parameters. In our case, the parameters are functions of e, so the list of 
parameters sometimes gets lengthy and the standard notation for these functions becomes 
cumbersome. We therefore introduce a more concise notation: 



Oil, 0(2, 



, Dtp 



z) = pF q (a 1 ,a 2 ,...,a p ;Pi,...,P q ;z) 



V A > • • • > Pq 

The generalized hypergeometric function has a power series representation: 
(oti,a 2 ,...,a p \ 

v V 

where (a)& is Pochhammer's symbol: 
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(C.126) 



(C.127) 



The power series converges for \z\ < 1. For z = 1, it converges if Res > 0, where 
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p-1 p 

i=l i=l 



(C.128) 



The hypergeometric function of type p+ \F q+ i has an integral representation in terms of the 
hypergeometric function of type p F q : 
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If a hypergeometric function has an upper and lower parameter that are equal, both param- 
eters can be deleted: 
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The simplest hypergeometric function is the one of type iFq. It can be expressed in an 
analytic form: 



iF (a; ;z) 



(1 



(C.131) 



The next simplest hypergeometric functions are those of type 2 F\. They satisfy transforma- 
tion formulas that allow an 2 F\ with argument z to be expressed in terms of an 2 F\ with 
argument z/(z — 1) or as a sum of two 2 F± 's with arguments 1 — z or 1/zor 1/(1 — z). The 
hypergeometric functions of type 2 Fi with argument z = 1 can be evaluated analytically in 
terms of gamma functions: 
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(C.132) 



The hypergeometric function of type ^F 2 with argument z = 1 can be expressed as a %F 2 
with argument z = 1 and different parameters 
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where s = Pi + (3 2 — ol\ — a 2 — a 3 . If all the parameters of a 3 F 2 are integers and half- 
odd-integers, this identity can be used to obtain equal numbers of half-odd-integers among 
the upper and lower parameters. If the parameters of a 3 F 2 reduce to integers and half- 
odd-integers in the limit e — > , the use of this identity simplifies the expansion of the 
hypergeometric functions in powers of e . 

The most important integration formulas involving 2 Fi hypergeometric functions is 
QU.129D with p = 2 and q = 1. Another useful integration formula is 
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(C.134) 
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This is derived by first inserting the integral representation for 2 iq in ( |C.129|) with integration 
variable t' and then evaluating the integral over t to get a 2-P1 with argument 1 + t'z. After 
using a transformation formula to change the argument to —t'z, the remaining integrals over 
t! are evaluated using ( |C.129|) to get 3-F 2 's with arguments —z. 

For the calculation of two-loop thermal integrals involving HTL averages, we require 
the expansion in powers of e for hypergeometric functions of type p F p ^i with argument 
1 and parameters that are linear in e. If the power series representation ( |C126| ) of the 
hypergeometric function is convergent at z — 1 for e = 0, this can be accomplished simply 
by expanding the summand in powers of e and then evaluating the sums. If the power series 
is divergent, we must make subtractions on the sum before expanding in powers of e. The 
convergence properties of the power series at z — 1 is determined by the variable s defined 
in ( |C128| ). If s > 0, the power series converges. If s — > in the limit e — > 0, only one 
subtraction is necessary to make the sum convergent: 
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If s — > — 1 in the limit e — > 0, two subtractions are necessary to make the sum convergent: 



cui, a 2 , . . . , a p 



Pi > • • • j A 



"p-i 



+ E 



T{(3 l )---T{(3^ l ) 
r(a 1 )r(a 2 ) • ■ - T{a vJ 
(ai) n (a 2 ) n ■ • ■ (a p ) n 



[((s + l)+tC(s + 2)} 



\ (^l)n ' ' ' C#?)« n! 

r(ai)r(a 2 ) • ■ -r(a p ) 



, (C.136) 



where t is given by 



= p (^-1X^-2) ^ (A-i)(A-2) 



(C.137) 
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The expansion of a P -F p _i hypergeometric function in powers of e is particularly simple if 
in the limit e — * all its parameters are integers or half-odd-integers, with equal numbers of 
half-odd-integers among the upper and lower parameters. If the power series representation 
for such a hypergeometric function is expanded in powers of e, the terms in the summand 
will be rational functions of n, possibly multiplied by factors of the polylogarithm function 
ip(n + a) or its derivatives. The terms in the sums can often be simplified by using the 
obvious identity 
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The sums over n of rational functions of n can be evaluated by applying the partial fraction 
decomposition and then using identities such as 
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The sums of polygamma functions of n + 1 or n + ^ divided by n + 1 orn+i can be evaluated 



using 
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where 71 is Stieltje's first gamma constant defined in flB.lOQ . The sums of polygamma 
functions of n + 1 or n + ~ can be evaluated using 
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We also need the expansions in e of some integrals of 2-F1 hypergeometric functions of y 
that have a factor of |1 — 2y\. For example, the following 2 integrals are needed to obtain 
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These integrals can be evaluated by expressing them in the form 



f Q dyy»-\l - yy-^l - 2y| F ^ 2 y^j = jf* ^^(l - 2 /)^ 1 (2 2/ - 1) F ^ 2 
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The evaluation of the first integral on the right side gives 3_F 2 hypergeometric functions with 
argument 1. The integrals from to | can be evaluated by expanding the power series 
representation ( |C.126| ) of the hypergeometric function in powers of e. The resulting series 
can be summed analytically and then the integral over y can be evaluated. 
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